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Abstract 

We have made next to leading order QCD fit to the deep inelastic 
spin asymmetries on nucleons and we have determined polarised quark 
and gluon densities. The functional form for such distributions was 
inspired by the unpolarised ones given by MRST (Martin, Roberts, 
Stirling and Thorne) fit. In addition to usually used data points (av- 
eraged over Q 2 ) we have also considered the sample containing points 
with the same x and different Q 2 . Our fits to both groups of data give 
very similar results. For the integrated quantities we get rather small 
gluon polarisation. For the non averaged data the best fit is obtained 
with vanishing gluon contribution at Q 2 = lGeV 2 . For comparison 
models with alternative assumptions about quark sea and in particular 
strange sea behaviour are discussed. 
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One has quite a lot of data on deep inelastic spin asymmetries for different 
nucleon targets. The data (recent and old) come from experiments made at 
SLAC [1-9], CERN [10-15] and DESY |T§, 0. The newest data on proton 



[§, [15], |17| and deuteron targets [§, || [L5[ have smaller statistical errors and 
hence dominate in \ 2 fits. 



Since the analysis of the EMC group results [[[(J an enormous interest 
were started in studying polarised structure functions. It was suggested Jl8 
that polarised gluons may be responsible for the little spin carried by quarks. 
Determination of polarised gluon distribution is particularly interesting in 
this context. After calculation of two loop polarised splitting functions |L9| 
several next to leading (NLO) order QCD fits were performed [20-24] and 
polarised parton distributions (i.e. for quarks, antiquarks and gluons) were 
determined. Many groups obtained rather highly polarised gluon contribu- 
tion (however this number is determined with a big error). The aim of this 
paper is to extend our next to leading order QCD analysis given in |25| by 
taking into account, in addition to all previously considered data, also pro- 
ton data from Hermes (DESY) and deuteron data from E155 experiment 



(SLAC). We will also use more recent MRST fit for partons [p6[] . In [25 



we got negligible gluon contribution at Q 2 = 1 GeV 2 . We will see that our 
result about gluon polarisation does not change very much. The method of 
choosing basic functions for fitting, used in this paper, differs from the other 
groups. As in |25j we will divide the data into two groups. Let us remind 



that many experimental groups published data ||, |15| sets for the near values 
of x and different Q 2 in addition to the " averaged" data where one averages 
over Q 2 (the errors are smaller and Q 2 dependence is smeared out). In prin- 
ciple when we take into account Q 2 evolution of polarised and unpolarised 
functions (in NLO analysis) the first group of data points, i.e. non averaged, 
should be considered. In most of the fits to experimental data only second 
group of data, namely with averaged Q 2 dependence, was used. We will 
make fits to the both sets of data (the first group contains 417 points and 
the second 159 points). The results for both groups of data are very similar 
(the same conclusion was already drawn in [ES ). We will also compare those 



results with the fits without Q 2 evolution taken into account (in other words 
assuming that asymmetries do not depend on Q 2 , as in our previous fits 
[p7[] ). One should add that many experimental groups have not succeeded in 
finding Q 2 dependence for approximately the same value of x and different 
Q 2 Wl O, |j"4"|. In our analysis we limit ourselves, as one usually does, to 
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the data with Q 2 > lGeV 2 . As was already mentioned in our earlier papers 
[ j27j] making a fit to spin asymmetries and not directly to gi(x, Q 2 ) enables 



to avoid the problem with higher twist contributions which are probably less 
important in such case ( see for example J28|) . Experiments on unpolarised 



targets provide information on the unpolarised quark densities q(x, Q 2 ) and 
G(x, Q 2 ) inside the nucleon. These densities can be expressed in term of 
q ± (x,Q 2 ) and G ± (x,Q 2 ), i.e. densities of quarks and gluons with helicity 
along or opposite to the helicity of the parent nucleon. The unpolarised 
quark densities are given by the sum of q + , q~ and G + , G~ , namely: 

q = q+ + q-, G = G + + G~ . (1) 

On the other hand the polarised DIS experiments give information about 
polarised parton densities, i.e. the difference of q + , q~ and G + , G~: 

Aq = q + -q~, AG = G + - G~. (2) 

We will try to determine q ± (x, Q 2 ) and G ± (x, Q 2 ), in other words, we will 
try to connect unpolarised and polarised data. In principle the asymptotic 
x behaviour of q + and q~ will be taken from the unpolarised data (up to the 
modifications when some leading order terms vanish). We will use the fit 
given by MRST p6| . In comparison with their previous fit (called R2) p9[ 
there is different behaviour at small x for quark and gluon distributions. It 
is of course very restrictive assumption that Aq and AG have the same be- 
haviour as q and G. On the other hand the small x behaviour of unpolarised 
structure functions is determined from the x values of Hera much smaller 
than in the polarised case. In our further analysis we will consider the in- 
tegrals over the region measured in the experiments with polarised particles 
with the hope that in this case the behaviour of q^ and G^ could be more 
plausible. The values of integrals in the whole region (0 < x < 1), involving 
asymptotic behaviour taken from the unpolarised structure functions, may 
be not as reliable as in the measured region. As an alternative we also use 
Regge type behaviour. 



It is known |20| that the behaviour of the quark and gluon distributions 
in small x region is extremely important in extrapolation of integrals over 
whole < x < 1 range. It could happen that in Aq = q + — q~ (when we 
assume that q + , q~ and q = q + + q~ have similar x dependence) most sin- 
gular x terms cancel (that is especially important in case of valence quarks 
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and sea contributions where the x behaviour is relatively singular). We will 
see later how such description infers the fits and calculated parameters. For 
the less singular distributions for Au v , Ad v and AM (total sea polarisation) 
there is no strong dependence of calculated quantities on the extrapolation 
in an unmeasured region but the fits have higher x 2 . One of the main tasks 
of considering NLO evolution in Q 2 is the determination of the gluon con- 
tribution AG. In MS scheme AG(x, Q 2 ) comes in through the higher order 
corrections. In our fits we obtain AG relatively small contrary to some ex- 
pectations. When we use the non-averaged sample of data the actually the 
best fit is with vanishing AG contribution. 

Let us start with the formulas for unpolarised quark parton distributions 
gotten (at Q 2 = 1 GeV 2 ) from the one of recent fits performed by Martin, 
Roberts, Stirling and Thorne |26j. For the valence quarks one get (in this fit 
one uses A^J 4 = 0.3 GeV and a s (Mf) = 0.120): 

u v (x) = 0.6051aT a5911 (i -x) 3 - 395 (l + 2.078^+ 14.56x), 

d v (x) = 0.0581x" a7118 (l-a;) 3 - 874 (l + 34.69v / ^ + 28.96x), (3) 

and for the antiquarks from the sea (the same distribution is for sea quarks) : 

2u(x) = 0.4M(x) -6(x), 

2d(x) = 0AM(x)+5{x), (4) 
2s(x) = 0.2M(x). 

In eq.(4) the singlet contribution M = 2[u + d + s] is: 

M(x) = 0.2004ar L2712 (l - x) 7 ' 808 (l + 2.283^ + 20.69x), (5) 

whereas the isovector part {5 — d — u) reads: 

8{x) = 1.290x 0183 (l - x) 9 - 808 (l + 9.987x - 33.34x 2 ). (6) 

For the unpolarised gluon distribution one gets: 

G(x) = 64.57aT a0829 (l - x) 6 ' 587 (l - 3.168v^ + 3.251x). (7) 

We assume, in an analogy to the unpolarised case, that the polarised 
quark distributions are of the form: x a (l — x) f3 P(y/x), where P{\/x) is a 



4 



polynomial in ^fx and the asymptotic behaviour for x — > and x — > 1 (i.e. 
the values of ct and /3) are the same (except for AM, see a discussion below) 
as in the unpolarised case. Our idea is to split the numerical constants 
(coefficients of P polynomial) in eqs.(3, 5, 6 and 7) in two parts in such a 
manner that the distributions q ± (x,Q 2 ) and G ± (x, Q 2 ) remain positive. At 
the end of the paper we will discuss the consequences of relaxing the positivity 
conditions. Our expressions for Aq(x) = q + (x)—q~(x) (q(x) = q + (x)+q~(x)) 
are: 

Au v (x) 
Ad v (x) 
AM(x) 
A5(x) 
AG(x) 

It is very important what assumptions one makes about the sea contri- 
bution. From the MRST fit for unpolarised structure functions the natural 
assumption would be: As = Ad/ 2 = Au/2. If we add the condition that 
SU(3) combination: as = Au + Ad — 2As should be equal to the value de- 
termined from the semileptonic hyperon decays, As is pushed into negative 
values and so is nonstrange sea. Instead of connecting As in some way to 
non-strange sea value we introduce additional free parameters for the strange 
sea contribution namely 

AM, = a-°- 7712 (l - xY- 808 (c ls + c 2sV ^). (9) 

For the strange quarks we have additional independent parameters. Hence, in 
our fits we will start with fourteen parameters. Comparing the expression (5) 
with (8) and (9) we see that in AM (and AM S ) there is no term behaving like 
2,-1.2712 a -j. sma }} x (hence, we assume that AM and hence all sea distributions 
have finite integral) which means that in AM coefficient in front of this term 
have to be splitted into equal parts in AM + and AM - (the most singular 
term in the sea contribution drops out). Hence, in the fitting procedure we 
are using functions that are suggested by the fit to unpolarised data. Maybe 
not all of them are important in the fit and it could happen that some of 
the coefficients in eqs.(8,9) taken as free parameters in the fit are small or in 
some sense superfluous. Putting them to zero or eliminating them increase 

5 



= x -°- 5911 (l - x y- 395 { ai + a 2 Vx~ + a 4 x), 

= x- a7118 (l - xf 87 \h + 6 2V ^ + b 3 x), 

= ^ a7712 (l-x) 7 - 808 (c 1 + c 2V ^), (8) 

= x ai83 (l - x) 9 - 808 c 3 (l + 9.987a; - 33.34x 2 ), 

= x -0.0829 (1 _ x) 6-587 {di + ^ + ^ 



X 2 only a little but makes x 2 /^df smaller. We will see that that is the 
case with some parameters introduced in eqs. (7,8). On the other hand we 
have still relatively strong singular behaviour of Au v and AM for small x 
values. For comparison we will also consider later the model in which most 
singular terms are put equal to zero namely = 0, which means 

that plus and minus components have the same coefficients for this power of 
x. In the less singular models the dependence of calculated parameters in 
the unobserved region (below x < 0.003) is weak. In the earlier papers we 
considered the extrapolation of various calculated integrals below x = 0.003 
up to assuming Regge type of behaviour for small x values. As will be 
discussed later the less singular models give significantly higher x 2 ■ 

In order to get the unknown parameters in the expressions for polarised 
quark and gluon distributions (eqs. (8,9)) we calculate the spin asymmetries 
starting from initial Q 2 = 1 GeV 2 for measured values of Q 2 and make a fit to 
the experimental data on spin asymmetries for proton, neutron and deuteron 
targets. The asymmetry Ai(x,Q 2 ) can be expressed via the polarised struc- 
ture function g\(x, Q 2 ) as 

^^.fi^sa.figfflp^^^,, (10) 

where R = [F 2 (l + 7 2 ) — 2xF 1 ]/2xF 1 whereas F\ and F 2 are the unpolarised 
structure functions and 7 = 2Mx/Q. We will take the new determined 
value of R from the [|30|. The factor (1 + 7 2 ) plays non negligible role for x 
and Q 2 values measured in SLAC experiments. Polarised structure function 
gi(x, Q 2 ) in the next to leading order QCD is related to the polarised quark 
and gluon distributions Aq(x, Q 2 ), AG(x,Q 2 ), in the following way: 

9i(x, Q 2 ) = \ J2 q e 2 g {Aq(x, Q 2 ) + ^[5c q * Aq(x, Q 2 ) 

+\5c g *AG(x,Q 2 )}} (11) 

with the convolution * defined by: 

(C*q)( Xl Q 2 )= f -C(-)q( Zl Q 2 ) (12) 

Jx Z Z 

The explicit form of the appropriate spin dependent Wilson coefficient Sc q 
and 5c g in the MS scheme can be found for example in ref. [|l]|. The NLO 
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expressions for the unpolarised (spin averaged) structure function is similar 
to the one in eq.(ll) with Aq(x, Q 2 ) — > q(x, Q 2 ) and the unpolarised Wilson 
coefficients are given in [[H], |32j . 

The Q 2 evolution of the parton densities is governed by the DGLAP 
equations |33]. For calculating the NLO evolution of the spin dependent par- 



ton distributions Aq(x, Q 2 ), AG(x, Q 2 ) and spin averaged q(x, Q 2 ), G(x, Q 2 ) 
ones we will follow the method described in |20], |32j . We will calculate Mellin 
n-th moment of parton distributions Aq(x, Q 2 ) and AG(x, Q 2 ) according to: 

Aq^ n \Q 2 ) = f 1 dxx^AqixiQ 2 ) (13) 
Jo 

and then use NLO solutions in Mellin n-moment space in order to calculate 
evolution in Q 2 for non-singlet and singlet parts. 

In calculating evolution of AS (n) (Q 2 ) and AG {n) (Q 2 ) with Q 2 we have 
mixing governed by the anomalous dimension 2x2 matrix . Having evolved 
moments one can insert them into the n-th moment of eq.(ll). 

9 { i\Q 2 ) = lE q e 2 q {Aq^(Q 2 ) + ^[5c^ ■ Aq^(Q 2 ) 

+I5c^-AG^(Q 2 )]} (14) 

and then numerically Mellin invert the whole expression. In this way we 
get g\(x,Q 2 ). The same procedure is applied for the unpolarised structure 
functions. Having calculated the asymmetries according to equation (10) for 
the measured in experiments value of Q 2 we can make a fit to a measured 
asymmetries on proton, neutron and deuteron targets. We will take into 
account all together 417 points (193 for proton, 171 for deuteron and 53 for 
neutron. We will use also the "experimental" value of as = 0.58 ±0.1 with 
enhanced (to 3a) error. 

The fit with all fourteen parameters from eqs.(8,9) gives \ 2 = 340.4. It 
seems that some of the parameters of the most singular terms are superfluous 
and we can eliminate them. We will put d\ = d,2 = (such assumption gives 
that 5G/G ~ x for small x), h\ = (the most singular term in Ad v ) and 
assume c\ s = c\ (i.e. the most singular terms for strange and nonstrange sea 
contributions are equal). Fixing these four parameters in the fit practically 
does not change \ 2 Du t improves \ 2 /Ndf- The resulting \ 2 P er degree of 
freedom is better than in the previous fit and one gets x 2 /Ndf — ifgiilo 
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=0.84. In this case we get the following values of parameters from the fit to 
all existing (above mentioned) data for Q 2 > lGeV 2 for spin asymmetries: 

a 4 = 15.7 ±0.42, 

(15) 



a i 


= 0.61 ±0.00, 


a 2 


= -6.1 ±0.19, 


b 2 


= -1.56 ±0.20, 


h 


= -0.43 ±0.49, 


Cl 


= -0.40 ±0.03, 


c 2 


= 4.15 ±0.00, 


Cls 


= c i, 


C2s 


= -0.28 ±0.83, 


c 3 


= -1.29 ±2.53, 






d 3 


= 2.01 ±11.2. 







Actually the parameter d 3 could be put equal to zero without increasing 
X 2 /N DF . We get in this case the smallest x 2 / n df =|§i^ = - 83 - That 
means that d% is not well determined in the fit and the best x 2 /^df is 
without gluonic contribution. 

The obtained quark and gluon distributions lead for Q 2 =1 GeV 2 to the 
following integrated (over x) quantities: Am = 0.80 ± 0.02, Ad = —0.65 ± 
0.03, As = -0.21 ± 0.05, Au v = 0.67 ± 0.02, Ad v = -0.59 ± 0.02, 2Au = 
0.14 ± 0.03, 2AJ = -0.07 ± 0.03. 

These numbers yield the following predictions: = Au ± Ad + As = 
-0.06 ± 0.07, a 3 = Au - Ad = 1.45 ± 0.02, AG = 0.04 ± 0.19, r? = 0.111 ± 
0.006, = -0.096 ± 0.006, Tf = 0.007 ± 0.005. 

We have positively polarised sea for up and negatively for down quarks 
and very strongly negatively polarised sea for strange quarks. Because of the 
big negative value of As the quantity ao is negative. The gluon polarisation 
is small. The value of 03 was not assumed as an input in the fit (as is the case 
in nearly all fits ||23|| ) and comes out slightly higher than the experimental 



value. A5, comes out relatively big from the fit (coefficient in front of A5 is 
equal to that in 5). 

The asymptotic behaviour at small x of our polarised quark distributions 
is determined by the unpolarised ones and these do not have the expected 
theoretically Regge type behaviour or pQCD which is also used by some 
groups, to extrapolate results to small values of x. Some of the quantities in 
our fit change rapidly for x < 0.003. 

Hence, we will present quantities integrated over the region from x=0.003 
to x=l (it is practically integration over the region which is covered by the 
experimental data, except of non controversial extrapolation for highest x). 
The corresponding quantities are Am = 0.85 (Au v = 0.56, 2Au = 0.29), 
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Ad = -0.48 (Ad v = -0.57, 2Ad = 0.09), As = -0.12, a = 0.25, AG = 
0.04, r? = 0.123, = -0.056, rf = 0.036, a 3 = 1.32. In this region the 
obtained values of sea contributions are relatively high and those of valence 
quarks relatively small. Gluon contribution practically vanishes. There is 
relatively strong dependence of different quantities in the unmeasured region 
(0 < x < 0.003). May be the unpolarised MRST parton distributions ( with 
the above mentioned modifications) do not describe quite correctly the small 
x behaviour of polarised parton distributions. On the other hand the fit to 
the data is very good. So, the values of integrated quantities in the measured 
region we consider as more reliable then in the whole region. With the value 
of As = —0.12 in the measured region of x we have Oq = 0.25 and with 
As = —0.21 in the whole region of x a Q becomes negative (-0.06). 

When we use the quantities calculated in the measured region and extend 
them to the full x region using asymptotic Regge behaviour for small x we 
get Am = 0.86 (Au v = 0.59, 2Au = 0.27), Ad = -0.51 (Ad v = -0.58, 
2Ad = 0.07), As = -0.14, a = 0.21, AG = 0.04, a 3 = 1.37. We have used 
x~ a behaviour for small x (with —0.25 < a < 0.25) and the quantities do not 
depend strongly on a specific value of a. For the given above values a = 
was used. 



HERMES E155 




0.05 0.1 0.2 0.5 0.02 0.05 0.1 0.2 0.5 



x x 

Figure 1: The comparison of our predictions for g± (x , Q 2 ) / Ff (x , Q 2 ) from 
basic fit with the recent averaged data for proton and deuteron targets. 

Now, we shall calculate T p , T n and T d in the measured region for Q 2 = 
5 GeV 2 and compare them with the quantities given by the experimental 
groups. We get in the region between x = 0.003 and x = 0.8 (covered by the 
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data) r? = 0.132 ± 0.006, T™ = -0.051 ± 0.007 and Tf = 0.037 ± 0.006. The 



experimental group SMC present [22] the following values in such region (for 
Q2 = 5 GeV 2 ); 

r; = o.i3o ±0.007, 

= -0.054 ±0.009, (16) 
Tf = 0.036 ±0.005. 

One can see that our results are in good agreement with experimental values. 
For comparison we have also made fits using formulas of the simple parton 
model (as in our papers before neglecting evolution of parton densities 
with Q 2 . More detailed result of these fits (integrated densities and so on) 
will be given later. 

In fig.l as an example we present our fit to the non averaged data in 
comparison with measured (averaged over Q 2 ) g±/Fi for new proton (Her- 
mes) and deuteron (E155) data. The curves are obtained by joining the 
calculated values of asymmetries corresponding to actual values of x and Q 2 
for measured data points. The curves are not fitted but the difference in 
fitted asymmetries for averaged and non-averaged data are very small. For 
asymmetries the curves with Q 2 evolution taken into account and evolution 
completely neglected do not differ very much so we do not present them. 

In figs. 2 and 3 we show the comparison of our predictions for g\ from the 
basic fit with the measured averaged values for proton, deuteron and neutron 
data. The values of g\ were calculated for the values of x and Q 2 measured for 
averaged data points in different experiments and then joined together. The 
agreement is good. On the other hand the spread of experimental points 
is still substantial. Polarised quark distributions for up and down valence 
quarks as well as non strange, strange quarks and gluons for Q 2 = 1 GeV 2 
are presented in figure 4. Dashed curves represent the ± and — components 
for different parton densities, the solid curves corespond to the difference 
of ± and — components, the sums of components (not shown) correspond 
to nonpolarised parton distributions. We see that especially polarised gluon 
distribution function is really tiny and does not resemble the distribution 
function for unpolarised case. 

This function is also quite different from the gluon distribution (given in 
0) used to estimate AG/G in COMPASS experiment planned at CERN 
p|. For x = 0.1 at Q 2 = 1 GeV 2 AG/G = 0.01 and is below a planned 
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Figure 2: The comparison of our predictions for gi(x, Q 2 ) with the measured 
structure functions in different experiments and on different nucleon targets. 
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Figure 3: The comparison of our predictions for g™(x, Q 2 ) calculated from 
basic fit with the measured structure functions in ElJf.2, E155 and Hermes 
experiments. 
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Figure 4: Our predictions for spin densities for quark and gluons at 
Q 2 = lGeV 2 . We present distributions for partons polarized along 
(xq + (x),xG + (x)) and opposite (xq~ (x) , xG~ (x) ) to the helicity of parent pro- 
ton. The polarized densities (i.e. the differences of above mentioned quanti- 
ties) are also shown. 
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experimental error. In fig. 5 we present spin densities for u quark. We show 
Au v (x, Q 2 ) obtained in the basic fit (solid line) as well as the same quantity 
from the fit with SU (3) symmetric sea (dashed curve) and from the fit where 
evolution in Q 2 was not taken into account (dotted line). In the second plot 
the same curves for Au(x, Q 2 ) are presented and one sees from that the values 
of Au from the basic fit and fit with SU (3) symmetric sea are very close in 
spite of the differences in the valence values. The similar conclusions are also 
true for the d quark densities. 




0.0050.01 0.05 0.1 0.5 1 0.0050.01 0.05 0.1 0.5 1 

x x 



Figure 5: Our predictions for xAu v (x,Q 2 = 1 GeV 2 ) and xAu(x,Q 2 = 
1 GeV 2 ) in different models. The solid curve is gotten using the parameters 
of our basic fit whereas dashed and dotted lines correspond to the quanti- 
ties calculated in the fit with a SU(3) symmetric sea and a fit without QCD 
evolution, respectively. 

Fixing the value of a§ is very important for the fit. When we relax the 
condition for a$ = 0.58 we get x 2 = 340.8, so this number practically does 
not change. We get the fit with the parameters not very different from our 
basic fit but with very small a 8 = 0.03 and positive As = 0.01. It seems that 
As is not well determined from the data on spin asymmetries alone but that 
does not influences strongly the values of Au, Ad and AC 

In order to check how the fit depends on the assumptions made about the 
sea contribution we have also made fit with Au = Ad = 2 As, the assumption 
that follows directly from MRST unpolarised fit. The x 2 value increases 
significantly and per degree of freedom one gets a number x 2 ' /N DF = 
=0.86 which is worse than in our basic fit. In this case we have Au = 0.80 
(Au v = 0.87, 2Au = -0.07), Ad = -0.61 (Ad v = -0.40, 2Ad = -0.21), 
As = —0.07, a = 0.11, AG = 0.07 and a 8 = 0.33. The quantity As must 
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be negative in order to get experimental value for ag and because of our 
assumption Au = Ad = 2As (the value of as does not come out correctly 
in the fit because of our assumption c ls = c\) we obtain negative values of 
non strange sea for up and down quarks. One sees that the values for sea 
polarisation depend very strongly on assumptions we made(in many papers 
|p0| , |2~T| , [34]] SU(3) symmetric sea is assumed that also together with fixing of 
a 8 value gives negative non strange sea). On the other hand Au = Au v + 2Au 
and Ad = Ad v + 2 Ad practically do not change (however, Au v and Ad v also 
change). Also AG does not change and is small. We get that the value of 
a§ is not coming out correctly from the fit to spin asymmetries. Fixing as 
and making specific assumption about As introduces shifts in nonstrange sea 
polarisation (and so in Au v and Ad v ) but Am, Ad do not change. Because 
the value of As is needed to determine a we decided to use additional free 
parameters for strange sea contribution in order to determine it (with fixed 
value of a 8 ) from the fit to experimental data. 

In many papers making fits to experimental data on spin asymmetries the 
assumption of SU(3) symmetric sea was made. We have also for comparison 
made fit with this assumption. The value of x 2 = 350.6 is substantially 
higher in comparison with our basic fit. In this case we have Am = 0.81 
(Au v = 0.87, 2Au = -0.06), Ad = -0.57 (Ad v = -0.40, 2Ad = -0.17), 
As = —0.11, do = 0.13, AG = 0.14, 03 = 1.38 and a§ = 0.47. In this case we 
also have shifts in values of valence and sea contributions similar to the case 
discussed above. 

Looking at the dependence of unpolarised quark and gluon densities we 
see that after elimination of most singular term in Ad v (x) the most singular 
behaviour for small x one has for Au v (x) and AM(x). For comparison we 
have investigated the model when in polarised densities these singular con- 
tributions are absent. In this case Au v and AM are ^Jx less singular than 
in our basic fit. For such a fit we get x 2 I^df = =0.87, i.e significantly 
higher than in our basic fit. We get in this case: Au = 0.77 (Au v = 0.57, 
2Au = 0.20), Ad = -0.38 (Ad v = -0.63, 2Ad = 0.25), As = -0.10, 
ao = 0.28, AG = 0.22. In such fit the integrated quantities taken over the 
whole range of < x < 1 and in the truncated one (0.003 < x < 1) differ 
very little. The quantity AG is positive and different from zero. So it is 
possible to get the fit with practically no change of integrated quantities in 
the region between x = and x = 0.003 but with significantly higher x 2 
value. For Q 2 =l GeV 2 we have r? = 0.122 and rj = -0.041. 
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The obtained results can be compared with the fit when instead of 417 
points for different x and Q 2 values we take spin asymmetries for only 160 
data points with the averaged Q 2 values for the same x (one has smaller 
errors in this case). In such fit the ratio of number of neutron to number 
of deuteron and proton data points is increased. It seems that the influence 
of neutron points is stronger than in basic fit x 2 I^df = j^zfo = 0-78 is a 
little bit better than in our basic fit. The integrated values for quark and 
gluon densities are:_ Aw = 0.79 {Au v = 0.65, 2Au = 0.14), Ad = -0.66 
(Ad v = -0.60, 2Ad = -0.06), As = -0.22, a = -0.09, AG = 0.31 and 
03 = 1.45. We see that averaging over Q 2 and different numbers of data 
points leads to very similar fit. The values for integrated valence densities 
and nonstrage sea contribution are only a bit shifted (Am and Ad in the whole 
region of integration do not differ from the basic fit for non averaged data and 
the same is also true for integrated quantities in the region 0.003 < x < 1). 
Integrated gluon density is relatively small and positive. A little bit higher 
value for AG = 0.31 ± 0.28 we do not consider as significant difference. 
X 2 /Ndf is very good and smaller then in [^4]] where the same experimental 
data sample and MRST parton distributions (modified for small values of x) 
were used. 

As was already mentioned before we have also made for comparison fits 
neglecting evolution of parton densities with Q 2 (formulas from the simple 
parton model). We get for non averaged data sample x 2 I^df = 4 3 1 4 8 9 l 9 9 =0.86 
(biger than in our basic fit: x 2 /Ndf =0.84): Am = 0.66 (Au v = 0.56, 
2Au = 0.10), Ad = -0.49 (Ad v = -0.49, 2Ad = 0.0), As = -0.20, a = 
-0.03, a 3 = 1.14, r? = 0.108, I™ = -0.082. For averaged data points we 
get x 2 /Ndf = 1 1 6 2 1 5 l 4 9 =0-83 (this number should be compared with x 2 I^df 
=0.78, the corresponding quantity from the NLO fit) and we have: Au = 
0.66 (Au v = 0.58, 2Au = 0.08), Ad = -0.48 (Ad v = -0.48, 2Ad = 0.0), 
As = —0.20, a = —0.03. Hence, x 2 P er degree of freedom is smaller in 
the case of averaged sample. We see that both fits give very similar results. 
It means that the averaging of data does not influence the fit when we do 
not take Q 2 evolution into account (the differences are also very small in the 
0.003 < x < 1 region). 

It has been pointed out [^l]] that the positivity conditions could be re- 
strictive and influence the contribution of polarised gluons. We have also 
made a fit to experimental data without such assumption for polarised par- 
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tons. The \ 2 value does not changed much x 2 I^df = 4ifzfo =0.84 and we 
get Am = 0.84 (Au v = 0.72, 2Au = 0.12), Ad = -0.74 (Ad v = -0.50, 
2Ad = -0.24), As = -0.24, a = -0.13, a 3 = 1.57, AG = 0.02. The results 
are a little bit different but the value of AG is not influenced by the positivity 
conditions. The same is also true in the case of averaged data. It seems that 
our positivity conditions are not very restrictive. 

We have made fits for two samples of data with averaged Q 2 values and 
non averged ones (adding neutron data from E154 and Hermes experiments) 
leading to very similar results for calculated parameters (except small differ- 
ence in AG). The value of a 3 was not fixed in the fit and comes out high 
in comparison with experimental value. In order to check the influence of 
different assumptions about strange sea we considered fits without fixing as 
value, with SU(3) symmetric sea and with modified sea contribution. The 
models with less singular behaviour for valence u quark and sea contribution 
were also discussed. In most of the modifications the x 2 value increases sig- 
nificantly. For comparison we have also considered fits to the simple parton 
model neglecting Q 2 dependence of parton densities. It seems that splitting 
of integrated densities Au, Ad into valence and sea contribution is model 
dependent (Au and Ad do not differ much). The integrated gluon contribu- 
tion comes out small. The best fits (measured by x 2 P er degree of freedom) 
we have for zero (for non averaged data points) or rather small (for averaged 
data) gluon polarisation. It seems that from results of our fits the perspective 
of measuring in COMPASS AG/G is not very encouraging. The experimen- 
tal accuracy still must be improved and probabely additional experiments 
are needed in order to make more precise statements about polarised quark 
and gluon densities. 
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